Physical condition and spin-resolved exchange correlation kernels 
in an inhomogeneous many electron system 
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We first exploit the spin symmetry relation fsS(0 — fis{~ C) f° r the exact exchange correlation 
kernel ffg(£) in an inhomogeneous many electron system with arbitrary spin polarization £. The 
physical condition required to satisfy the specific symmetry relation ffg(Q = /fs C (C) is derived and 
examined for simple ferromagnetic-nonmagnetic structure by taking the electrochemical potential 
into account. The condition is then applied to several composite systems useful in spintronics 
applications such as the magnetic system with net spin polarization. 
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Recently the spin transport properties have been at- 
tracting great interest for their potential applications to 
spintronics and quantum computation [lj, |2j, y, U| ■ One 
of the central problems in these fields is controlling the 
distribution of spin-polarized carriers in multicomponent 
structure alternating ferromagnetic and nonmagetic ma- 
terials. In this system, the exchange correlation (XC) 
kernel can be useful for giving direct insight into the car- 
rier population in spin-polarized channels. Exchange cor- 
relation kernel (XCK) /* s c , (£) is denned by 



_ d 2 E xc (r,r';() 
dn s ' (r)dn s (r') 



(1) 



where E xc is the total XC energy functional in many elec- 
tron system with spin polarization £. /* S C /(C) is a basic 
concept in describing many body correlation effects in an 
inhomogeneous electron liquid and satisfies the symme- 
try relation f£g(() = C) where s(s) refers to the 
majority (minority) spin. The spin symmetry relation 
of XCK plays a significant role to understand the na- 
ture of the spin-spin response function (which contains 
spin symmetric and anti-symmetric parts of XCK) in ar- 
tificial composite structure. Unfortunately, E xc , a key 
ingradient of XCK, is not capable of providing accurate 
ffg, in spite of numerous studies including density gra- 
dient corrections [3, |a> 13 • They fail to reduce the mean 
absolute error to the desired level in the chemical bonding- 
energies. That means that the desired chemical accuracy 
has not been reached yet. While the "mixed scheme" 
HI El Ei EH combining the density functional theory 
(DFT) with other methods such as quantum Monte Carlo 
simulations and coupled cluster calculations applied sep- 
arately to the short and long range parts of the electron- 
electron interaction has been proposed as an alternative, 
the short range part of the Coulomb interaction between 
electrons is still well described by semilocal functionals 
0,0]. On the other hand, gradient corrected density 



functionals (GCDF) HIE0 

have been used for studies of 
electronic structures but mainly been restricted to unpo- 
larized systems. Calculation of density functionals in the 
spin-polarized system has been extended to the local spin 
density functional scheme combining with other approx- 
imation methods |lC|. Extension of GCDF to the spin- 
polarized system is not available yet. Hence the sym- 
metry relation of the "exact" spin-resolved XCK has not 
been well established in an inhomogeneous spin polarized 
system. In this paper, we first exploit the symmetry rela- 
tion of the "exact" spin-resolved XCK and the condition 
required to satisfy the specific symmetry relation (SSR) 
fS(C) = /£(0 of the exact XCK. 

Although considerable researches have been devoted to 
investigate spin current I s tg), there has been no attempt 
to interpret XCK by directly measurable quantities such 
as I s (s) °t — 3/ g^ s) where /i s (s) is the electrochemical po- 



dx 



tential (ECP) in spintronics. In view of the fact that 
the required condition to satisfy SSR could be related 
to the spin density variation Vn s (j) and then Vn s (|j to 
the ECP variation V/i s ( g ), information on a system can 
be obtained straight by experimental observations of spin 
related phenomena. Also, theoretical ECP can exactly be 
checked through the symmetry relation of XCK. Hence 
we first propose the proper situations satisfying SSR in 
alternating multilayer system. We also give the accurate 
relation of corresponding spin-resolved pair correlation 
functions g 3s '(r) and sw(r'). 

GCDF is given by the sum of the kinetic energy T of a 
nonintcracting particle system, potential energy U, and 
unknown functional E xc 

E[n(r)]=T[n{r)} + U[n(r)}+E* c , (2) 

where one particle density n(r) is written by 



n(r) = N 



*0isi, • • • , r N s N )\ 2 dr 2l • • • , drjv(3) 



2 



with spins of N electrons. 

In pair density theory giving more accurate value of 
the ground state energy than one particle density [l3| . 
the spin-summed pair density is given by n(r!,r 2 ) = 
N(N-i) J2 SiS2 7 SlS2 (ri,r 2 ). Here lsiS2 (n, r 2 ) is the 
spin-resolved diagonal of the two-body reduced density 
matrix |l2l | 



ls 1 s 2 {r 1 ,r 2 ) 



, rArSAr)| 2 dr 3 , 



N- 



(4) 



The exact energy density functional is given by E[n(x)] = 
T[n(x)] + U with x = (ri, r 2 ). Here, E/ consists of the ex- 
ternal potential of a given pair and the interaction poten- 
tial between particles forming two pairs at xt — (r^i, r.; 2 ) 
and Xj = (rji,rj2) 0|. Hence, the total XC energy 
functional E* c in interacting system can be described in 
terms of pair density n(x) 

E xc [n(r)} = T[n{x)] - T[n(r)] + U[n(x)] - U[n(r)]. (5) 

E KC contains the correlated kinetic term T xc = T[n(x)] — 
T[n(r)] as well as the interparticle interaction potential. 
Here, T xc denotes the difference between the exact ki- 
netic term T[n(x)] on interacting scheme and noninter- 
acting counter part T[n(r)]. / x |, the second derivatives 
of E KC with respect to spin densities n s (r) and n s (r') at 
a given pair position Xi — (r, r'), is now written by 

pcfn = VrVr-[r xc [n(x)] + U[n(x)} - U[n(r)}} 
Js ~ s ^ Vn s (r)Vn s (r') ' [ 1 

For the exactly defined XCK, the symmetry relation 
fss(() = fss(~0 is trivial in an inhomogeneous broken 
spin symmetry system since majority and minority spins 
interchange their orientations and positions with the re- 
versed polarization — £. That is, based on the assumption 
that the physical condition 



Vn s (r / ) _ Vn 5 (r) 
Vn,(r') ~ Vn s (r) 



(7) 



is fulfilled (i.e., the ratios of spin density gradients are 
the same at two different positions r and r'), SSR 



d 2 E xc {r,r';() _ d 2 E KC (r , r' '; C) 
dn s (r)dn s (r') dn s (r)dn,s(r') 



f*:(o (8) 



is obtained trivially. In other words, only ratios of spin 
density variations at different sites are required to investi- 
gate SSR. This condition is valid in various density vary- 
ing systems. By examining the condition given by Eq.JTJ 
in various spin valve systems, we can investigate the va- 
lidity of SSR. Given the relation between ECP fi s /s) and 
the nonequilibrium spin carrier density n s (j) in metal- 
lic and nonmetallic system, spin density variation Vn s (s) 
can be obtained. 



In a highly degenerate system, the density variations 
for spin-up and spin-down carriers are gjyen, in the pres- 
ence of an electric field E — — V$, by 



Vn 



( 5 ) = eVL> s(g) (e F )[M s (s) +e$] + eD s(g) (e F )[V/i s(s) 

where D s tg\(ep) is the spin- up (spin-down) density of 
states at the Fermi level. Depending on the dimension of 
multilayer structure, D s r s \ is varied but the gradient of 
D s /g)(ep) at the fixed Fermi level vanishes always. Hence 
the condition J2J can be written by 



eE] 



D s (e F ) Vn,-eE 



D s (e F ) V/i« - eE 



D s {e F ) Vn s -eE 



D s (e F ) V/i« - eE 



(9) 



In nonmetallic region of a homogcmeous system with no 
space charge, the ratios of local variances Vn s /Vn 5 are 
the same. In doped systems, spin polarization can be 
created keeping the total number of electrons and holes 
constant, 8n s + Sn s — For a constant equilibrium 
spin density n , it is a trivial situation satisfying SSR 

Vn°+V5n s ' 



-1. 



We first consider the general ECP /i s (g) 14] in alter- 
nating ferromagnets (F) and nonmagnets (N) to see in 
what region SSR is satisfied. In a homogeneous F with 
"up" magnetization, the local variation of fj, s (g) with re- 
spect to x perpendicular to the layer is derived, at zero 
temperature, by 

- cE , (P+V rrXn)^ 



dx 



(10) 



where (3 and l^* are the bulk spin asymmetry coefficient 
(— 1 < /3 < 1) and spin diffusion length, respectively, in 

(n) 

F. The constants K- are determined from the proper 
boundary conditions in the nth layer. For a given bulk 
spin resistivity p s r s -\ — l/c s (s) = 2p F [l =p P] in F, can 
be determined from the spin-scattering measurements 
[Til Hi|. Here, p F is the fixed resistivity in F obtained 
by measurement and — (+) sign in front of (3 corresponds 
to the up (down) spin. Near the interface, scattering is 
localized significantly in an interfacial region. Since the 
spin flip scattering rate is not easily manipulated and 
identified near the interface, obtaining accurate f3 is not 
easy. As the stabilization of (3 is determined by the im- 
purity concentration, [3 should be checked with care at 
positions r and r' in a region especially near the inter- 
face. From Eqs.(JZJ) and 0, the condition required to 
satisfy SSR is 



(3+1 



(3+1 



(3-1 



(»-e->r 

hf 



l sf 



)■ (11) 



For a homogeneous N with "up" magnetization, the 
local variation of the general ECP in the nth layer is 
derived, similarly, by 



1 



dx 



= eE± w [K 2 



(12) 
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The ratio of spin density variations is simply — and 
SSR is trivially satisfied. 

Based on the results of the general cases given by 
Ecis. H10|) and (|12(l . we consider a simple F/N structure 
as schematically shown in Fig. 1(a). A spin-polarized cur- 
rent with density j s ^ flows from F (x < 0) into N (x > 0) 
along the direction perpendicular to the interface. p s (s) 
can be expressed in terms of the spin-resolved conduc- 
tivity <7 S ( S ) and current density j s ( S ) with the spin accu- 
mulation balanced by the spin flip scattering. The local 
variation of ECP is written by 9l " " — 



in F and 



dx 



2e Pi<sjs(s) m N where is the re- 



sistivity of N [16J. In this case, the physical condition, 
Eq.Ql becomes 



(1 - 0)j s + TE 



(1 + (3)j s + TE 



(1 - 0)j a + TE 



in F and 



3. + 



3s + 



(1 + 0)j 3 + TE 



3s + TE 



h + ^e 



(13) 



(14) 



in N with T = 1/(2/9*), (i—F, N). Far from the interface 
(x 3> 0) in N with equilibrium spin state, cr s = erg and 
j s = j s flfij ]. This is an obvious situation satisfying SSR 
since Vn s /Vn s is constant in the region x ^> 0. Near the 
interface with /i s — fi s ^ (due to the spin scattering), 
the conditions l|13f) and l|14(l should be checked with care 
in order to satisfy SSR. 

For the case of F/N/F structure as shown in Fig. 1(b), 
ECP is an odd function of x [l7L Hfij . In this system, 
P state or AP state exists since the spin orientations in 
two F's are either parallel(P) or antiparallel(AP). In each 
state, ECP is different and we can obtain the required 
condition by considering each case. In AP state, the 
ECP gradients are derived by 



+ ^cosh^)^)], -f < x < 



( dx 

I aft 



x<-f 



(15) 

where J is the total current density and op and o>j are the 
total conductivities in F and N. In P state, parameters 
Ai, A2, Bi and cosh(-rf-) of Eq. l|15|l are substituted by 

A[, —A' 2 , B[ and sinh(™-). In the region x < — the 

...*/. 
condition required to satisfy SSR is 



°JL 

0". 



(16) 



In the region — ^ < x < 0, it becomes simply = 
— 1. In two dimensional electron gas, the conductivity 
is proportional to the density of states at the constant 
Fermi level so that the conductivity cr s ( S ) in N between 
two ferromagnets satisfies a s ( S ) = 0"n/2 where <tn is the 
total conductivity in N. While V/z s 7^ Vfi s in P state, 
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FIG. 1: Schematic illustration of spin valve system, (a) a 
single ferromagnetic layer (F) and nonmagnetic layer (N) (b) 
F/N/F structure for parallel spin arrangement between fer- 
romanetic layers. A polarized spin current is injected into 
nonmagnet via ferromanet perpendicular to the layer, (c) 
F1/N1/F2/N2 structure consisting of thin enough Nl and F2 
and infinitely thick Fl and N2. (d) F/N1/N2 structure with 
a highly doped semiconductor Nl between F and N2. 



V/i s = V/ig in AP state in N. That results an unpolarized 
current flow in AP state and a polarized current flow 
(is ^ js) in P state. 

Let us take a limiting case of 1^ 3> ijj such that 
T s f = 00 in N. In this limit, one can approximate 
fx s (s) — Hos(s) + r^s)^ [13 j expanding up to the first 
order in x. The first term /j,q s (s) denotes ECP with no 
spin effect. The constant value of means that SSR 



is also satisfied in this limit. Hence the required condi- 
tion is trival to satisfy the symmetry relation of XCK 
given by Eq. (JHJ ■ The assumption of r s / = 00 is valid in 
several materials such as Si-doped GaAs Q • Experimen- 
tal observations in spin valve systems composed of these 
materials can be analyzed newly according to SSR. 

The dilute magnetic semiconductor (DMS) has been 
suggested as an alternative to F due to large conduc- 
tivity mismatches between F and N. In DMS system, 
the spin polarization results from the indirect electron- 
electron interaction. The ECPs of DMS/N/DMS struc- 
ture 0| satisfy Ea. lji"5|) . except that parameters A4, B i: 
<7f, and Iff are to be replaced respectively by Ai, Bi, 
ctd, and iff which are determined by the proper bound- 



ary conditions. Hence the investigation of SSR in this 
system is similar to that of F/N/F structure. 

In F1/N1/F2/N2 multilayer as shown in Fig. 1(c), elec- 
trodes are attached to Fl and N2 to measure the voltage 
V across the multilayer. In the ith layer, ECP is given 
by H3 



±x/hsf 



(17) 



where the parameters Ci and A is ^ are functions of (3 and 
have different values in different parts of the multilayer. 
From the general solution of ECP, we see that ^± = 
in F and ^7 = —1 in N. When the ratios ^— at r and 
r' are the same, SSR is always satisfied. 

Let us consider a metallic junction [^], in which ECPs 
have the same form as in F1/N1/F2/N2 structure. For 
the case of Py/Cu/Py system, A is ^ S ) and k s f can be ob- 
tained in experiments. For example, when polarized car- 
riers are injected into the Cu wire via the Py pad, Aj s (g) 
are given by T-A*i(l i ap y )/2 for each part i of the struc- 
ture with Ipysf = Ap y = 2nm and ap y = 0.2 Then, 
^7 = —1.5 at different sites in the present case. 

We consider an F/N1/N2 structure, shown in Fig. 1(d), 
composed of F (x < 0), Nl (0 < x < xq) with the conduc- 
tivity <7ni, and N2 (x > xq) with <7N2- A highly doped 
semiconductor Nl is often placed between N2 and F. The 
spin densities are given by Q 



±[A e 



-x/ld 



(x~x )/l u 



?[1 — 9ss(r)] H2]. Hence, Eq.Q can be written by 



n s Vg ss (r) + n s Vg sS (r) 



n s Vg S s{r) + n s Vg Ss (r) 



n s Vg ss (r) + n s Vg sS {r)} 



n s Vg ss {r) + n s Vg Ss (r) 



(18) 



(20) 

The spin-resolved pair correlation functions including the 
correlation effect can not be evaluated accurately using 
quantum Monte Carlo algorithm 23]. That is, the accu- 
rate condition given by Ea. Q20[l can be used as an indi- 
cator testing the accuracy of the spin-resolved pair cor- 
relation functions. 

In summary, we have investigated the symmetry re- 
lation of the "exact" spin-resolved XCK in broken spin 
symmetry system. We have shown the proper cases sat- 
isfying SSR in multicomponcnt structure. Only by using 
ECP gradient proportional to the measured spin current, 
the properties of XCK can be easily checked. The proper 
condition 22°. I _ SJM. I can ^ e a standard of the cor- 

V n <rjr V™<r W 

rectness of spin-related measurements especially in the 
trivial situations, for example, a homogeneous or non- 
magnetic system. Hence, we have proposed new method 
to interpret not only spin current but also the properties 
of XCK directly in spintronics. We also give the accu- 
rate relation of spin-resolved pair distribution functions 
which can also be used to test the precision of the pair 
distribution functions. 
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for < x < xq and 



s(s) 



±[A 2 e 



-(x-x )/l d ] 



(19) 



for x > xo where l u and Id are up- and down-stream spin 
diffusion lengths. When Vn°| r = 0, it is a trivial situation 
for SSR to be satisfied in undoped and doped systems. 

The SSR for spin currents flowing through the quan- 
tum dots can also be examined in a similar way. In 
magnetic dots, it is interesting to examine the in- 
direct exchange interaction between conduction elec- 
trons and localized spins due to magnetic impurities. 
In such magnetic system involving localized d elec- 
trons, the spin density is written by n s /g\ (r) = 2Ja qp 
f m 2 ^EiF(2k F (\r - Ri\) < Sf > Q where n c , 
Ep, Ri, and ft^ are, respectively, the total number of 
conduction electrons, unperturbed Fermi energy, sites 
of magnetic ions, and the volume of a unit cell. Here, 
_ zcosz-sim; ^ Q n tmg accoun ^ ) S y S t e m belongs 

to an inhomogeneous spin polarized system and SSR is 
satisfied for a certainty from the fact that ^7 = — 1. 

For all these situations satisfying SSR, the properties of 
the pair correlation function can also be known. Spin de- 
pendent eletron density in a homogeneous spin polarized 
system is given, in general, by n s (r) = n s [l — g ss (r)\ + 
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